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I. INTRODUCTION 



The long-sought Higgs boson (s) h of electroweak symmetry breaking in particle physics 
may soon be observed at the CERN Large Hadron Collider (LHC) through the diphoton 
decay mode (h — > 77). Purely hadronic standard model processes are a copious source 
of diphotons, and a narrow Higgs boson signal at relatively low masses will appear as a 
small peak above this considerable background. A precise theoretical understanding of the 
kinematic distributions for diphoton production in the standard model could provide valuable 
guidance in the search for the Higgs boson signal and assist in the important measurement 
of Higgs boson coupling strengths. 

In this paper we address the theoretical calculation of the invariant mass, transverse mo- 
mentum, rapidity, and angular distributions of continuum diphoton production in proton- 
antiproton and proton-proton interactions at hadron collider energies. We compute all con- 
tributions to diphoton production from parton-parton subprocesses through next-to-leading 
order (NLO) in perturbative quantum chromodynamics (QCD). These higher-order contri- 
butions are large at the LHC, and their inclusion is mandatory for quantitatively trust- 
worthy predictions. We resum initial-state soft and collinear logarithmic terms associated 
with gluon radiation to all orders in the strong coupling strength a s . This resummation 
is essential for physically meaningful predictions of the transverse momentum (Qt) distri- 
bution of the diphotons at small and intermediate values of Qt, where the cross section is 
large. In addition, we analyze the final-state collinearly-enhanced contributions, also known 
as 'fragmentation' contributions, in which one or both photons are radiated from final-state 
partonic constituents. We compare the results of our calculations with data on isolated 
diphoton production from the Fermilab Tevatron The good agreement we obtain with 
the Tevatron data adds confidence to our predictions at the energy of the LHC. The present 
work expands on our recent abbreviated report 0], and it may be read in conjunction with 
our detailed treatment of the contributions from the gluon-gluon subprocess [3||. 

Our attention is focused on the production of isolated photons, i.e., high-energy photons 
observed at some distance from appreciable hadronic remnants in the particle detector. 
The rare isolated photons tend to originate directly in hard QCD scattering, in contrast to 
copiously produced non-isolated photons that arise from nonperturbative processes such as 
7r and r] decays, or from via quasi-collinear radiation off final-state quarks and gluons. 

We evaluate contributions to continuum diphoton production from the basic short- 
distance channels for 77 production initiated by quark-antiquark and (anti) quark- gluon 
scattering, as well as by gluon-gluon and gluon- (anti) quark scattering proceeding through 
a fermion-loop diagram. At lowest order in QCD, a photon pair is produced from qq anni- 
hilation [Fig. CD (a)]. Representative next-to-leading order (NLO) contributions to qq + qg 
scattering are shown in Fig.[T](b)-(e). They are of 0(a s ) in the strong coupling strength 0,0]. 
Production of 77 pairs via a box diagram in gg scattering [Fig. Q] (h)] is suppressed by two 
powers of a s compared to the lowest-order qq contribution, but it is enhanced by a product 
of two large gluon parton distribution functions (PDFs) if typical momentum fractions x 
are small [6fl. The 0(a^ s ) or NLO corrections to gg scattering include one-loop gg — > 77$ 
diagrams (i) and (j) derived in Ref. @, S|, as well as 4-leg two-loop diagrams (1) computed 
in Ref. 0, H3|. In this study we also include subleading contributions from the process (k), 
9(ls ~^ HQs via the quark loop, where qs = J2i =u d s (<?i + <?«) denotes the flavor-singlet 
combination of quark scattering channels. 
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Figure 1: Representative partonic subprocesses that contribute to continuum diphoton production. 
All leading-order and next-to-leading order direct production subprocesses, i.e., contributions (a)- 
(e) and (h)-(l), are included in this study. Diagrams (f) and (g) are examples of single-photon one- 
and two-fragmentation. 



Factorization is a central principle of hadronic calculations in perturbative QCD, in which 
a high-energy scattering cross section is expressed as a convolution of a perturbative partonic 
cross section with nonperturbative parton distribution functions (PDFs), thus separating 
short-distance from long-distance physics. The common factorization is a longitudinal no- 
tion, in the sense that the convolution is an integral over longitudinal momentum fractions, 
even if some partons in the hard-scattering process have transverse momenta that border the 
nonperturbative regime. Unphysical features may then arise in the transverse momentum 
(Qt) distribution of a color-neutral object with high invariant mass (Q), such as a pair of 
photons produced in hadron-hadron collisions. When calculated in the common factorization 
approach at any finite order in perturbation theory, this distribution diverges as Qt — ► 0, 
signaling that infrared singularities associated with Qt —> have not been properly iso- 
lated and regulated. These singularities are associated with soft and collinear radiation from 
initial-state partons shown by the diagrams in Figs. CD (b), (d), and (i). 

A generalized factorization approach that correctly describes the small-Qr region was 
developed by Collins, Soper, and Sterman (CSS) [H| and applied to photon pair produc- 
tion [z,[H, 13]. In this approach the hadronic cross section is expressed as an integral over 
the transverse coordinate (impact parameter). The integrable singular functions present 
in the finite-order differential distribution as Qt — > are resummed, to all orders in the 
strong coupling a s , into a Sudakov exponent, and a well-behaved cross section is obtained 
for all Qt values. As explained in Sec. II, our resummed calculation is accurate to next-to- 
next-to-leading-logarithmic (NNLL) order. It is applicable for values of diphoton transverse 
momentum that are less than the diphoton mass, i.e., for Qt < Q. When Qt ~ Q, terms 
of the form ln n (Q T / Q) become small. A perturbative expansion with a single hard scale is 
then applicable, and the cross section can be obtained from finite-order perturbation theory. 
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In addition to the initial-state logarithmic singularities, there is a set of important final- 
state singularities which arise in the matrix elements when at least one photon's momentum 
is collinear to the momentum of a final-state parton. They are sometimes referred to as 
'fragmentation' singularities. At lowest order in a s , the final-state singularity appears only 
in the qg — > 77^ diagrams, as in Fig. [U (e). There are various methods used in the literature 
to deal with the final-state singularity, including the introduction of explicit fragmentation 
functions D 1 (z) for hard photon production, where z is the light-cone fraction of the interme- 
diate parton's momentum carried by the photon. These single-photon "one-fragmentation" 
and "two-fragmentation" contributions, corresponding to one or both photons produced in 
independent fragmentation processes, are illustrated by the diagrams in Figs. [Q (f ) and (g). 
In addition, a fragmentation contribution of entirely different nature arises when the 77 pair 
is relatively light and produced from fragmentation of one parton, as discussed in Sees. Ill C 21 
and lHI A~3l A full and consistent treatment of the final-state logarithms beyond lowest order 
would require a joint resummation of the initial- and final-state logarithmic singularities. 

In the work reported here, we are guided by our interest in describing the cross section 
for isolated photons, in which the fragmentation contributions are largely suppressed. A 
typical isolation condition requires the hadronic activity to be minimal (e.g., comparable to 
the underlying event) in the immediate neighborhood of each candidate photon. Candidate 
photons can be rejected by energy deposit nearby in the hadronic calorimeter or the presence 
of hadronic tracks near the photons. A theory calculation may approximate the experimental 
isolation by requiring the full energy of the hadronic remnants to be less than a threshold 
"isolation energy" Ef° in a cone of size AR around each photon. The two photons must be 
also separated in the plane of the rapidity 77 and azimuthal angle <p by an amount exceeding 
the resolution Ai? 77 of the detector. The values of Ef°, AR, and Ai? 77 serve as crude 
characteristics of the actual measurement. The magnitude of the final-state fragmentation 
contribution depends on the assumed values of Ef°, AR, and AR J7 . 

An additional complication arises when the fragmentation radiation is assumed to be 
exactly collinear to the photon's momentum, as implied by the photon fragmentation func- 
tions Dj(z). The collinear approximation constrains from below the values of z accessible 
to D 1 (z): z > z min . The size of the fragmentation contribution may depend strongly on the 
values of Elf and z min as a consequence of rapid variation of D^(z) with z. 

In our work we treat the final-state singularity using a prescription that reproduces desir- 
able features of the isolated cross sections while bypassing some of the technical difficulties 
alluded to above. For Qt > Ef°, we avoid the final-state collinear singularity in the qg 
scattering channel by applying quasi-experimental isolation. When Qt < Elf. , we apply 
an auxiliary regulator which approximates on average the full NLO rate from direct qg and 
fragmentation cross sections in this Qt range. Two prescriptions for the auxiliary regulator 
(subtraction and smooth-cone isolation inside the photon's isolation cone) are considered 
and lead to similar predictions at the Tevatron and the LHC. 

We begin with our notation in Sec. Ill Al followed by an overview of the procedure for 
resummation of initial-state multiple parton radiation in Sec. Ill Bl The issue of the final- 
state fragmentation singularity is discussed in Sec. Ill CI Our approach is compared with that 
of the DIPHOX calculation [lj], in which explicit fragmentation function contributions are 
included at NLO, but all-orders resummation is not performed. Our theoretical framework 
is summarized in Sec. Ill Dl 

In Sec. II I II we compare the predictions of our resummation calculation with Tevatron 



4 



data. Resummation is shown to be important for the successful description of physical Qt 
distributions, as well as for stable estimates of the effects of experimental acceptance on 
distributions in the diphoton invariant mass. We compare our results with the DIPHOX 
calculation [iJ] and demonstrate that the requirement Qt < Q further suppresses the effects 
of the final-state fragmentation contribution, beyond the reduction associated with isolation. 
Next, we present our predictions for distributions of diphoton pairs produced at the energy 
of the LHC. Various distributions of the diphoton pairs produced from the decay of a Higgs 
boson are contrasted with those produced from QCD continuum processes at the LHC, 
showing that enhanced sensitivity to the signal can be obtained with judicious event selection. 
Our conclusions are presented in Sec. 



II. THEORY OVERVIEW 
A. Notation 

We consider the scattering process h\(P\) + /^(Pa) — ► l{Pz) + l{P&) + X, where h\ and 
/12 are the initial-state hadrons. In terms of the center-of-mass collision energy \/S, the 
invariant mass Q, transverse momentum Qt, and rapidity y of the 77 pair, the laboratory 
frame momenta P± and P^ of the initial hadrons and g M = P% + Pf of the 77 pair are 

^{1,0,0,1}; (1) 
{1,0,0,-1}; (2) 

(3) 




1 Q 2 + Q\ cosh y,Q T ,0, \J Q 2 + Q T sinh y j 



The light-cone momentum fractions for the boosted 2 — > 2 scattering system are 



_ 2(P 2|1 ■ q) _ y/Q T +me* v m 



Decay of the 77 pairs is described in the hadronic Collins-Soper frame [15|. The Collins- 
Soper frame is a rest frame of the 77 pair (with g M = {Q, 0,0,0} in this frame), chosen 
so that (a) the momenta Pi and P2 of the initial hadrons lie in the Oxz plane (with zero 
azimuthal angle), and (b) the z axis bisects the angle between Pi and — P%. The photon 
momenta are antiparallel in the Collins-Soper frame: 

P£ = — {0, sin 9* cos 99* , sin sin 99* , cos^*} , (5) 

Pf = ^ {0, — sin 6* cos — sin 6* sin (p*, — cos 8*} , (6) 

where 6* and tp* are the photon's polar and azimuthal angles. In this section, we derive 
resummed predictions for the fully differential 77 cross section da/(dQ 2 dydQ T dQ*), where 
dfl* = dcosOxdcp* is a solid angle element around the direction of P 3 in the Collins-Soper 
frame defined in Eq. (JSJ). The angles in the Collins-Soper frame are denoted by a "*" 
subscript, in contrast to angles in the lab frame, which do not have such a subscript. The 
parton momenta and helicities are denoted by lowercase pi and Aj, respectively. 
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B. Resummation of the initial-state QCD radiation 



For completeness, we present an overview of the finite-order and resummed contributions 
associated with the direct production of diphotons. At the lowest order in the strong coupling 
strength a s , photon pairs are produced with zero transverse momentum Qt- The Born 
qq — ► 77 cross section corresponding to Fig. [1] (a) is 



d<7qq 



dQ 2 dy dQ^dfl* 



KQt) S ' a* U/hi ( X l , f*F)f§i/ha fa , Hf) , (7) 



Born i=u,u,d,d,- 



where f qi /h(x, Hf) denotes the parton distribution function (PDF) for a quark of a flavor i, 
evaluated at a factorization scale //p of order Q. The prefactor 

W*M - (0) 1 + C0S2 6 * /oN 

with 

jo) = a 2 {Q)etn , , 

1 ~ 2N C Q* ' { ' 

is composed of the running electromagnetic coupling strength a = e 2 /47r evaluated at the 
scale Q, fractional quark charge e, = 2/3 or —1/3, and number of QCD colors N c = 3. 

The lowest-order gg — > 77 scattering proceeds through an amplitude with a virtual quark 
loop (a box diagram) shown in Fig. [TJ (h) . Its cross section takes the form 



da gg 



= S(Q T ) ^ 9 fg/fo (Zl , Vf) fg/ha (X 2 ,Hf), (10) 
Born ^ 



dQ 2 dy dQ T dQ* 
where the prefactor 

E g (Q = afL g {6*) (11) 

depends on the polar angle 6>* through a function L g (9*) presented explicitly in Ref. [3]. 
The overall normalization coefficient 



(0) a 2 (Q)a 2 s (Q) 
9 " 327rQ 2 (iV c 2 - 




(12) 



involves the sum of the squared charges e 2 of the quarks circulating in the loop. 

The NLO direct contributions, represented by Figs, [p (b)-(e), (i)-(l) and denoted as 
P(Q,Q T ,y,^*), are computed in Refs. @, |, H, 0, M, @, E9|. The NLO 2^3 differen- 
tial cross section grows logarithmically if the final-state parton is soft or collinear to the 
initial-state quark or gluon, i.e., when Qt of the 77 pair is much smaller than Q. These 
"initial-state" logarithmic contributions are summed to all orders later in this subsection. 
The NLO qg cross section also contains a large logarithm when one of the photons is pro- 
duced from a collinear q ) — >V7 splitting in the final state. This "final-state" collinear limit 
is discussed in Section Hi CI 



6 



With contributions from the initial-state soft or collinear radiation included, the NLO 
cross section is approximated in the small-Qr asymptotic limit by 



A q q(Q, Qt, V, fi*) 



E 



S 



6(Q T )F i>5 (Q,y,9*)+F ii+ (Q,y,Q T ) 



(13) 



i=u,u,d,d,.. 

in the qq + qg scattering channel, and by 

1 



Agg(Q,Q T ,y,^*) = ^ S 5 



5(Q T )F giB {Q,y,9.) + F gt+ (Q,y,Q T ) 
Z , g (9.,<p.)F' g {Q,y,Q T ) 



(14) 



in the gg+gqs scattering channel. The functions F a ^{Q, y, 9*) and F^ + (Q, y, Qt) for relevant 
parton flavors a are listed in Appendix O They include 'plus function' contributions of 
the type [Q^ 2 h"i p (Q 2 /Qt)] + with p > 0, universal functions describing soft and collinear 
scattering, and process-dependent corrections from NLO virtual diagrams. 

The qq + qg asymptotic cross section A qq -(Q, Qt, y, fi*) is proportional to the angular 
function £j(#*), the same as in the Born qq — > 77 cross section, cf. Eq. (JJJ). Similarly, the 
gg + gqs asymptotic cross section A gg (Q, Q T , y, includes a term proportional to the Born 
angular function Yi g (9*). In addition, A gg (Q,QT,y,Q*) contains another term proportional 
to Eg(#*, v 9 *) = L' g (9*) cos 2^*, where L' g (9*) is derived in Ref. jjjj. This term arises due to the 
interference of Born amplitudes with incoming gluons of opposite polarizations and affects 
the azimuthal angle (97*) distribution of the photons in the Collins-Soper frame (§]. 

The small-Qr representations in Eqs. ( 1X31 ) and ( 1X41 ) can be used to compute fixed-order 
particle distributions in the phase-space slicing method. In this method, we choose a small 
Qt value Q^ p in the range of validity of Eqs. (fX3|) and (fl4"J) . If the actual Qt in the 
computation exceeds Qt P , we calculate the differential cross section using the full 2^3 
matrix element. When Qt is smaller than Qt P , we calculate the event rate using the small- 
Qt asymptotic approximation A(Q, Q T , y, f)*) and 2^2 phase space. Hence, the lowest 
bin of the Qt distribution is approximated in the NLO prediction by its average value in 
the interval < Qt < Qt P , computed by integration of the asymptotic approximations. 

The phase-space slicing procedure is sufficient for predictions of observables inclusive in 
Qt, but not of the shape of da/dQT distributions. The latter goal is met by all-orders summa- 
tion of singular asymptotic contributions with the help of the Collins-Soper-Sterman (CSS) 
method H, 0, 13]. The small-Qy resummed cross section is denoted as W(Q, Qt, y, 0*) 
and given by a two-dimensional Fourier transform of a function W(Q, b, y, fi*) that depends 
on the impact parameter b: 



W(Q,Q T ,y,n* 



db 

db 
(2<p 



e^ rb W(Q,b,y,Q,) 

e l ^ rb w pert {Q, K y, nje-ww. 



(15) 



In this equation, W(Q, b, y, £7*) is written as a product of the perturbative part 
W pert (Q, K,y, Q*) and the nonperturbative exponent exp (—Fnp(Q, b)) , which describe the 
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dynamics at small (6 < 1 GeV 1 ) and large (b > 1 GeV impact parameters, respectively. 
The purpose of the variable 6* is reviewed below. 

If Q is large, the perturbative form factor W pert dominates the integral in Eq. (TT5l) . It is 
computed at small b as 

w pert {QAv,e*) = E ^^(^*) e ~ 5a(Q,6) 

a 

X [C a/ai ® / ai / fcl ] Ol, b] n) [Ca/ a2 ® /aa/fcj (^2, 6; A*)- (16) 

The "hard-vertex" function E a (#*)/i 2 ((5, 0*) is the normalized cross section for the Born scat- 
tering aa — > 77, with a — u,u, d, d, ... in qq — > 77, and a = a = g in gg — > 77. The Sudakov 
exponent 



«S a (Q,6) 



Cf/6 2 ^ 2 



A (Ci, fi) In ( ) + B a (d, C 2 , /i) 



(17) 



is an integral of two functions A a (Ci,p,) and B a (Ci,C 2 ,/2) between momentum scales C\/b 
and C 2 Q, and Ci and C 2 are constants of order c = 2e _7B = 1.123... and 1, respectively. The 
symbol [C a / ai ® / tt i/h] (a?, b; /i) stands for a convolution of the kx— integrated PDF f ai /h(%, A 4 ) 
and Wilson coefficient function C / 0l (x, 6; C1/C2, /i), evaluated at a factorization scale /x and 
summed over intermediate parton flavors ai: 



[Ca/ai ®/ai/fc] A*) = ^2 J -J C a/ax [j^^^J ^i/^' 

0-1 



A 4 ) 



(18) 



We compute the functions fr a , *4. a , i3 a and C a / ai up to orders a s , , a 2 , and a s , respectively, 
corresponding to the NNLL accuracy of resummation. The perturbative coefficients at these 
orders in a s are listed in Appendix lAl 

The subleading contribution from the nonperturbative region b > 1 GeV" 1 is included in 
our calculation using a revised "6*" model [18J, which provides excellent agreement with px- 
dependent data on Drell-Yan pair and Z boson production. In this model, the perturbative 
form factor W pert (Q, &*, y, f2*) in Eq. (fl5l) is evaluated as a function of fo* = bj (1 + 6 2 /b^^) 1 / 2 , 
with b m ax = 1.5 GeV" 1 . The factorization scale \i in [C <S> f] is set equal to c$\/b~ 2 + Q 2 ni , 
where Qj„i is the initial scale of order 1 GeV in the parameterization employed for f a /h(x, aO, 
for instance, 1.3 GeV for the CTEQ6 PDFs Q. We have W pert (K) = W pert (b) at b 2 < 6 2 



J max ' 



and Wpert(&*) = Wpertibmax) at fe 2 3> Hence, this ansatz preserves the exact form of the 

perturbative form factor W per t(Q,b,y,Q,*) in the perturbative region of small 6, while also 
incorporating the leading nonperturbative contributions (described by a phenomenological 
function J^NpiQ^b)) at large b. 

The form of Tnp{Q, b) found in the global pt fit in Ref. [l8| suggests approximate inde- 
pendence of J^NpiQ, b) from the type of qq scattering process. It is used here to describe the 
nonperturbative terms in the leading qq — »■ 77 channel. We neglect possible corrections to 
the nonperturbative contributions arising from the final-state soft radiation in the qq chan- 
nel and additional \fS dependence affecting Drell-Yan-like processes at x < 10~ 2 |20(] , as 
these exceed the accuracy of the present measurements at the Tevatron. The experimentally 
unknown J-"np(Q, b) in the gg channel is approximated by !Fnp(Q, b) for the qq channel, 



multiplied by the ratio Ca/Cf = 9/4. This choice is motivated by the fact that the lead- 
ing Sudakov color factors in the gg and qq channels are proportional to Ca = 3 and 
Cf = 4/3, respectively. The uncertainties in the 77 cross sections associated with Fnp(Q, b) 
are investigated numerically in Ref. 0|. 

In the region Qt ~ Q, collinear QCD factorization at a finite fixed order in a s is ap- 
plicable. In order to include non-singular contributions important in this region, we add to 
W(Q, Qt, y, ^*) the regular piece Y(Q, Qt, y, £1*), defined as the difference between the NLO 
cross section P(Q,QT,y,Q*) and its small-Qr asymptotic approximation A(Q, Qt, y, ^*): 

^QdQ^Tydn = W (Q>QT,y,n*) + P(Q,Q T ,y,^) - A(Q,Q T ,y,n*) 

= W(Q, Q T , y, 0.) + Y(Q, Q T , y, O.)- (19) 

At small Qt, subtraction of A(Q, Qt, y, fi*) in Eq. (fT9l) cancels large initial-state radia- 
tive corrections in P(Q,Qt,V,Q*), which are incorporated in their resummed form within 
W(Q, Qt, V, 0*) ■ At Qt comparable to Q, A(Q,Q T ,y,Q*) cancels the leading terms in 
W(Q, Qt, y, but higher-order contributions remain from the infinite tower of logarith- 
mic terms that are resummed in W. In this situation the W + Y cross section drops below 
the finite-order result P(Q, Qt, V, ^*) at some value of Qt (referred to as the crossing point) 
in both the qq + qg and gg + gqs channels, for each Q and y. We use the W + Y cross section 
as our final prediction at Qt values below the crossing point, and the NLO cross section P 
at Qt values above the crossing point. 

A few comments are in order about our resummation calculation. The hard- vertex contri- 
bution Ti a (9^)h^(Q, 6*) and the functions B a {C\, C 2 , p) and C a / ai (x, b; C/C 2 ,fi) can be varied 
in a mutually compensating way while preserving the same value of the form factor W up to 
higher-order corrections in a s . This ambiguity, or dependence on the chosen "resummation 



scheme" [2l| within the CSS formalism, can be employed to explore the sensitivity of the- 
oretical predictions to further next-to-next-to-next-to-leading logarithmic (NNNLL) effects 
that are not accounted for explicitly. 

The perturbative coefficients in Appendix[A]are presented in the CSS resummation scheme 



111 ], our default choice in numerical calculations, and in an alternative scheme by Catani, 
de Florian and Grazzini (CFG) (HJ. In the original CSS resummation scheme, the B and 
C functions contain the finite virtual NLO corrections to the 2^2 scattering process, 
whereas in the CFG scheme the universal B and C depend only on the type of incident 
partons, and the process-dependent virtual correction is included in the function h a . The 
difference between the CSS and CFG schemes is numerically small in 77 production at both 
the Tevatron and the LHC [!]. 

In the gg + gqs scattering channel, the unpolarized resummed cross section includes an 
additional contribution from elements of &/r-dependent PDF spin matrices with opposite 
helicities of outgoing gluons [!]. The NLO expansion of this spin-flip resummed cross section 
generates the term proportional to £'(#*,</?*) oc cos2</?* in the small-Q^ asymptotic cross 
section, cf. Eq. (TT4l) . Although the logarithmic spin-flip contribution must be resummed in 
principle to all orders to predict the dependence in the gg + gqs channel, it is neglected 
in the present work in view of its small effect on the full 77 cross section. 

When integrated over Qt from to scales of order Q, the resummed cross section be- 
comes approximately equal to the finite-order (NLO) cross section, augmented typically by 
a few-percent correction from integrated higher-order terms logarithmic in Qt- Inclusive 
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observables that allow such integration (e.g., the large-Q region of the 77 invariant mass 
distribution) are approximated well both by the resummed and NLO calculations. How- 
ever, the experimental acceptance constrains the range of the integration over Qt in parts 
of phase space and may break delicate cancellations between integrable singularities present 
in the finite-order differential distribution. In this situation (e.g., in the vicinity of the kine- 
matic cutoff in dcr/dQ discussed in Sec. UTTl l the NLO cross section becomes unstable, while 
the resummed cross section (free of discontinuities) continues to depend smoothly on kine- 
matic constraints. We see that the resummation is essential not only for the prediction of 
physical Qt distributions in 77 production, but also for credible estimates of the effects of 
experimental acceptance on distributions in the diphoton invariant mass and other variables. 

C. Final-state photon fragmentation 

1. Single-photon fragmentation 

In addition to the QCD singularities associated with initial-state radiation [described by 
the asymptotic terms in Eqs. (fT3l) and (fill)], other singularities arise in the 0{a s ) process 
<l(.Pi) + 9{P2) — > 7(^3) +7(^4) +q(j>5) [Fig. ED (e)] when a photon is collinear to the final-state 
quark. In this limit, the qg — > ^77 squared matrix element grows as l/s 7 5, when s y5 — > 0, 
where s 7 5 is the squared invariant mass of the collinear jq pair. In this limit, the squared 
matrix element factors as 

\M(qg -> <m)| 2 « ^P^ q (z)\M(qg -> gi )\ 2 (20) 

into the product of the squared matrix element \M.{qg — > q-f) \ 2 for the production of a 
photon and an intermediate quark, and a splitting function P 7 ^ q (z) = (1 + (1 — z) 2 )/z for 
fragmentation of the intermediate quark into a collinear 75 pair. In Eq. (l20l) z is the light- 
cone fraction of the intermediate quark's momentum carried by the fragmentation photon, 
a nd etj is the charge of t he intermediate quark. When the photon-quark separation Ar = 
\l (jfc ~ V-y) 2 + (<P~5 ~ ¥i) 2 m the plane of pseudorapidity rj = — log(tan(6'/2)) and azimuthal 
angle cp in the lab frame is small, as in the collinear limit, s 7 5 ~ E^E^Ar 2 , where Et^ 
and Et5 are the transverse energies of the photon and quark, with Et = -Esin^. Note that 
Et5 — Qt at the order in a s at which we are working. Therefore, contributions from the 
final-state collinear, or fragmentation, region are most pronounced at small Ar and relatively 
small Qt- 1 

A fully consistent treatment of the initial- and final-state singularities would require a 
joint initial- and final-state resummation. In the approaches taken to date, the fragmentation 
singularity may be subtracted from the direct cross section and replaced by a single-photon 

"one- fragmentation" contribution q + g — > (q 7) +7, where "(g — > 7)" denotes collinear 
production of one hard photon from a quark, described by a function D^z^) at a light- 
cone momentum fraction z and factorization scale /1. Single-photon "two-fragmentation" 



1 In the soft, or E5 — > 0, limit, the final-state collinear contribution is suppressed, reflecting the absence of 
the soft singularity in the qg — > 977 cross section. 
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contributions arise in processes like g + g — > (q 7) + (q ^—^ 7) and involve convolutions 
with two functions D 1 (z,fi) (one per photon). The lowest-order Feynman diagrams for the 
one- and two-fragmentation contributions are shown in Figs. QJf) and QJg), respectively. 
Parameterizations must be adopted for the nonperturbative functions D^(z, u) at an initial 
scale fi — jiQ. This is the approach followed in the DIPHOX calculation [14j], in which the 
sum of real and virtual NLO corrections to direct and single-7 fragmentation cross sections 
is included. When explicit fragmentation function contributions are included, the inclusive 
rate is increased by higher-order contributions from photon production within hadronic jets. 
However, much of the enhancement is suppressed by isolation constraints imposed on the 
inclusive photon cross sections before the comparison with data. Nevertheless, fragmentation 
contributions surviving isolation may be moderately important in parts of phase space. 
An infrared-safe procedure can be formulated to apply isolation cuts at each order of 



a 



22|, |23|, |24|] . This procedure encounters difficulties in reproducing the effects of isolation 
on fragmentation contributions, because theoretical models reflect only basic features of the 
experimental isolation and may introduce new logarithmic singularities near the edges of the 
isolation cones. 

As mentioned in the Introduction, the magnitude of the fragmentation contribution de- 
pends on the values of isolation parameters E™°, AR, and Ai? 77 , modeled only approximately 
in a theoretical calculation. The collinear approximation constrains from below the values 
of z accessible to D 7 (z,iS): z > z min = (1 + E'^ / E Tl )~ x . If D y (z, (J,) varies rapidly with z, 
the fragmentation cross section is particularly sensitive to the assumed values of Elf° and 
z min . For instance, if D 7 (z,/i) ~ 1/z, the fragmentation cross section is roughly proportional 
to Et° under a typical condition E™ / 'Er 1 <C 1. Such nearly linear dependence on E^° of 
the fragmentation cross section da/dQr is indeed observed in the DIPHOX calculation, as 
reviewed in Sec. IHIl In reality, some spread of the parton radiation in the direction trans- 
verse to the photon's motion is expected. The treatment of kinematics in parton showering 
programs like PYTHIA results in somewhat different dependence on z [12| compared to the 
collinear approximation, hence in a different magnitude of the fragmentation cross section. 

In this work we adopt a procedure that reproduces desirable features of the isolated cross 
sections, while bypassing some of the difficulties summarized above. To simulate experimen- 
tal isolation, we reject an event if (a) the separation Ar between the final-state parton and 
one of the photons is less than AR, and (b) Ets of the parton is larger than E^°. This 
condition is applied to the NLO cross section P(Q,Qt,u,£1*), but not to W(Q,QT,y,Q*) 
and A(Q, Qt, y, #*), as these correspond to initial-state QCD radiation and are free of the 
final-state collinear singularity. 

This quasi-experimental isolation excludes the singular final-state direct contributions at 
E T5 > E™ and Ar < AR (or s 75 < E T ^E Tb AR 2 ). It is effective for Q T > E™°, but 
the collinear direct contributions survive when Qt < E™°. The integrated (but not the 
differential) fragmentation rate in the region Qt < Elf° may be estimated from a calculation 
with explicit fragmentation functions. In our approach, we do not introduce fragmentation 
functions, but we apply an auxiliary regulator to the direct qg cross section at Qt < E^° 
and Ar < AR. In our numerical study we find that this prescription preserves a continuous 
differential distribution except for a small finite discontinuity at Qt = E^°. It approximately 
reproduces the integrated qg rate obtained in the DIPHOX calculation at small Q Tl for the 
nominal E^°. 
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Figure 2: Lowest-order Feynman diagrams describing fragmentation of the final-state partons into 
photon pairs with relatively small mass Q. 



Two forms of the auxiliary regulator are considered below, based on subtraction of the 
leading collinear contribution and smooth-cone isolation [25fl. In the first case, we subtract 
the leading part Eq. (1201) of the direct qg matrix element when Et5 < E^° and Ar < AR. 
We take z — 1 — p s -Ps/ips 'Pf +Ps • Pb + Pf -Pb), where p^, p$, and p^ are the four-momenta 



of the fragmentation photon, fragmentation quark, and spectator photon, respectively [26]. 
This prescription is used in most of the numerical results in this paper. 

In the second case, we suppress fragmentation contributions at Ar < AR and Et5 < E^° 
by rejecting events in the AR cone that satisfy Et5 < x(Ar), where \(Ar) is a smooth 
function satisfying x(0) = 0, x(A-R) = E™ . This "smooth-cone isolation" [25|| transforms 
the fragmentation singularity associated with D^(z, /jl) into an integrable singularity, which 
depends on the assumed functional form of x(Ar). The cross section for direct contributions 
is rendered finite by this prescription without explicit introduction of fragmentation functions 
D 7 (z,fi). For our smooth function, we choose x(Ar) = E^°(l— cos Ar) 2 /(1— cos AR) 2 , which 
differs from the specific form considered in Ref. [25|], but still satisfies the condition x(0) = 0. 
Our earlier results in Ref. [2] are computed with this prescription. Here we employ it only 
in a few instances for comparison with the subtraction method and obtain similar results. 

Differences between the two prescriptions can be used to quantify sensitivity of the pre- 
dictions to the treatment of the Qt < E^° and Ar < AR region. The two prescriptions 
yield identical predictions outside of this restricted region, notably at Qt > E™ , where 
our NLO perturbative expression P(Q, Qt, U, ^*) in the qq + qg channel is controlled only 
by quasi-experimental isolation and coincides with the corresponding direct cross section in 
DIPHOX. The default subtraction prescription predicts a vanishing da/dQT in the extreme 
Qt — > limit, while the smooth-cone prescription has an integrable singularity in this limit, 
avoided by an explicit small-Qr cutoff in the calculation of our Y-piece. Both prescriptions 
are free of the logarithmic singularity at Qt = E^° arising in the fixed-order (DIPHOX) 
calculation. 



2. Low-Q diphoton fragmentation 

Another class of large radiative corrections arises when the 77 invariant mass Q is smaller 
than the 77 transverse momentum Qt- In this case, one final-state quark or gluon fragments 

into a low-mass 77 pair, e.g. as q + g — * (q ^—i 77) + g. The lowest-order contributions 
of this kind are shown in Fig. El The process is described by a 77-fragmentation func- 
tion D 11 (zi, Z2, n), different from the single-photon fragmentation function D 7 (z, fi). This 



12 



new "two-photons from one-fragmentation" contribution is not included yet in existing cal- 
culations, even though similar fragmentation mechanisms have been studied in large-Qr 
Drell-Yan pair production 27, H^]. The importance of \ow-Q 77-fragmentation may be el- 
evated in some kinematic regions for typical experimental cuts. They can be removed by 
adjustments in the experimental cuts, as discussed in Sec. | 



D. Summary of the calculation 

We conclude this section by summarizing the main features of our calculation. Full direct 
NLO cross sections, represented by the graphs (a)-(e), (h)-(l) in Fig. [TJ are computed, and 
their initial-state soft/collinear logarithmic singularities are resummed at small Qt in both 
the qq + qg and gg + gqs channels. The perturbative Sudakov functions A and B and Wilson 
coefficient functions C in the resummed cross section W are computed up to orders af, a 2 s , 
and a s , respectively, corresponding to resummation at NNLL accuracy. 

Our resummation calculation requires an integration over all values of impact parameter 
b, including the nonperturbative region of large b. In our default calculation of the resummed 



cross section, we adopt the nonperturbative functions introduced in Ref. [I8j. We consider 
two resummation schemes, the traditional scheme introduced in the CSS paper as well as an 
alternative scheme [HI]. The comparison allows us to estimate the magnitude of yet higher- 
order corrections that are not included. The size of these effects is different in the qq + qg 
and gg + gqs channels but not particularly significant in either (§]. 

The final-state collinear singularity in the qg scattering channel is avoided by applying 
quasi-experimental isolation when Qt > EJf° and an auxiliary regulator when Qt < E™ to 
approximate on average the full NLO rate from direct qg and fragmentation cross sections 
in this Qt range. Two prescriptions for the auxiliary regulator (subtraction and smooth 
isolation inside the photon's isolation cone) are considered and lead to similar predictions at 
the Tevatron and LHC. 

The singular logarithmic contributions associated with initial-state radiation are sub- 
tracted from the NLO cross section P to form a regular piece Y, which is added to the 
small-Qr resummed cross section W to predict the production rate for small and intermedi- 
ate values of Qt- In the gg + gqs channel, we also subtract from P a new singular spin-flip 
contribution that affects azimuthal angle dependence in the Collins-Soper reference 
frame. We switch our prediction to the fixed-order perturbative result P at the point in Qt 
where the cross section W + Y drops below P. This crossing point is located at Qt of order 
Q in both qq + qg and gg + gqs channels. 



III. COMPARISONS WITH DATA AND PREDICTIONS 

Our calculation of the differential cross section da / (dQdQTdydQ*) is especially pertinent 
for the transverse momentum Qt distribution in the region Qt < Q-, for fixed values of 
diphoton mass Q (cf . Section IIII A II) . It would be best to compare our multiple differential 
distribution with experiment, but published collider data tend to be presented in the form 
of singly differential distributions in Q, Qt, and Aip = </? 3 — </?4 in the lab frame, after 
integration over the other independent kinematic variables. We follow suit in order to make 
comparisons with Tevatron collider data, but we recommend that more differential studies 
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be made, and we comment on the features that can be explored. We show results at the 
energy of the Tevatron collider and then make predictions for the Large Hadron Collider. 

The analytical results of Sec. [Til are implemented in our computer code. As a first step, 
resummed and NLO 77 cross sections are computed on a grid of discrete values of Q, Qt, 
and y by using the resummation program Legacy described in Refs. 29|, 3(3] ■ At the second 



stage, matching of the resummed and NLO cross sections is performed, and fully differential 
cross sections are evaluated by Monte-Carlo integration of the matched grids in the latest 



version of the program ResBos [31|, ]32j. The calculation is done for Nf = 5 active quark 
flavors and the following values of the electroweak and strong interaction parameters [331 ]: 

G F = 1.16639 x 10" 5 GeV" 2 , m z = 91.1882 GeV, (21) 
a(m z ) = 1/128.937, a s (m z ) = 0.1187. (22) 

The following choices of the factorization constants are used: C\ = C 3 = 2e _7B ~ 1.123..., 
and C2 = C4 = 1. The choice C4 = 1 implies that we equate the renormalization and 
factorization scales to the invariant mass of the photon pair, = fip = Q, in the fixed- 
order and asymptotic contributions P(Q,Qt,V,Q*) and A{Q, Qt, y, fi*). We use two-loop 
expressions for the running electromagnetic and str ong couplings a(/i) and «s(/i), as well as 
the NLO parton distribution function set CTEQ6M [19J with Q ini = 1.3 GeV. For calculations 
with explicit final-state fragmentation functions included, we use set 1 of the NLO photon 
fragmentation functions from Ref. [33 ]. 



A. Results for Run 2 at the Tevatron 



1. Kinematic constraints 

In this section, we present our results for the Tevatron pp collider operating at y/~S = 
1.96 TeV. In order to compare with the data from the Collider Detector at Fermilab (CDF) 
collaboration [l|, we make the same restrictions on the final-state photons as those used in 
the experimental measurement (unless stated otherwise): 

transverse momentum p T > p Tmin = 14 (13) GeV for the harder (softer) photon, (23) 
and rapidity |?/ 7 | < 0.9 for each photon. (24) 

We impose isolation conditions described in Section III CI assuming the nominal isolation 
energy E™ = 1 GeV specified in the CDF publication, along with AR = 0.4, and Ai? 77 = 
0.3. 

We also show predictions for the constraints that approximate event selection conditions 



used by the Fermilab D0 Collaboration [35] : p T > p Tmin = 21 (20) GeV for the harder 
(softer) photon, |y 7 | < 1.1, and E^°/E T = 0.07 for each photon, for the same AR and Ai? 77 
values as in the CDF case. 

A scatter plot of event distributions from our theoretical simulation for CDF kinematic 
cuts and arbitrary luminosity is shown in Fig. [3l The events are plotted versus the invariant 
mass Q, transverse momentum Qt, rapidity separation \Ay\ = \yhard ~ Vsoftl, an d azimuthal 
separation Aip = \<fhard — Vsoftl (with < Aip < it) between the harder and softer photon 
in the lab frame, as well as the cosine of the polar angle 9* in the Collins-Soper frame. It 
can be seen from the figure that Aip is correlated with the difference Qt — Q- Events with 
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Kinematical distributions in the Tevatron Run-2 (CDF) 
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Figure 3: The diphoton event distribution from the theoretical simulation for y/S = 1.96 GeV, with 
the selection criteria imposed in the CDF measurement, as a function of the various kinematic 
variables described in the text, shown for Qt < Q and Qt > Q separately. 

Qt < Q (Qt > Q) tend to populate regions with A<p > n/2 (Atp < 7r/2). The extreme case 
Qt = relevant to the Born approximation corresponds to A(p = tt. 

The p T cuts suppress the mass region Q < 2 \J VrminVrmin ~ 27 GeV at Atp w n and Qt < 
25 GeV at Atp w 0, leading to the appearance of a kinematic cutoff in the invariant mass 
distribution and a "shoulder" in the transverse momentum distribution, as shown in later 
sections. Our theoretical framework is applicable in the region Qt < Q (large A<p), where the 
dominant fraction of events occurs. The appearance of singularities in the NLO calculation 
at Qt — * and the fact that there are two different hard scales, Qt and Q, relevant for the 
event distributions in the 1ow-Qt region require that we address and resum large logarithmic 
terms of the form \og(Q /Qt). Different and interesting physics becomes important in the 
complementary region Q T > Q (small Atp), a topic we address in Sec. IIII A3l 

2. Tevatron cross sections 

We compare our resummed and finite-order predictions for the invariant mass (Q) dis- 
tribution of photon pairs, shown in Fig. 0] as solid and dashed lines, respectively. The 
finite-order cross section is evaluated at 0(a s ) accuracy in the qq + qg channel and at 0(af) 
accuracy in the gg + gqs channel. These finite-order calculations are performed with the 
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Figure 4: Invariant mass distributions of photon pairs in pp — > 77 JT at y/S = 1.96 TeV with 
QCD contributions calculated in the soft-gluon resummation formalism (red solid) and at NLO 
(blue dashed). The calculations include the cuts used by the CDF collaboration whose data are 
shown 

phase-space slicing method described in Sec. Ill Bl When integrated over all Qt, as in the 
da/dQ distribution at large Q, the resummed logarithmic terms from higher orders in a s 
produce a relatively small NNLO correction, such that the resummed and finite-order mass 
distributions in Fig. [Hare close to one another in normalization and shape. Both distributions 
also agree with the CDF data in this Q range within experimental uncertainties. 

The shape of da/dQ at small Q is affected by the cuts in Eq. ( 1231 ) on the transverse 
momenta p T of the two photons. In addition to being responsible for the characteristic 
cutoff at Q ~ 27 GeV explained in the previous subsection, the cuts on the individual 
transverse momenta p T also introduce a dependence of the invariant mass distribution on 
the shape of the Qt spectrum of the 77 pairs. Because of this correlation between the Q 
and Qt distributions, the discontinuities in da/dQT as Qt — » 0, when computed at finite 
order, make finite-order predictions for da/dQ somewhat unstable. 

The finite-order expectation for the transverse momentum distribution da/dQT (i-e., the 
integral of P(Q, Qt, V, ^*) over Q, y, and or P for brevity) is shown as a dashed curve 
in Fig. E](a). It exhibits an integrable singularity in the small-Qr limit. Terms with inverse 
power and logarithmic dependence on Qt, associated with initial-state radiation as Qt — >■ 0, 
are extracted from P and form the asymptotic contribution, denoted as A (dotted curve). 
In the figure, both P and A are truncated at a small value of Qt, that is, not drawn all 
the way to Qt = 0. The curves for P and A are close at small values of Qt, signaling that 
the initial-state logarithmic singularities dominate the NLO distribution. The difference Y 
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Figure 5: Transverse momentum distributions in pp — ► 77X at y/S = 1.96 TeV along with the CDF 
data: (a) the fixed-order prediction P (dashes) and its asymptotic approximation A (dots); (b) the 
full resummed cross section (solid), obtained by matching the resummed W + Y to the fixed-order 
prediction P (dashed, same as in (a)) at large Qt- 



between the P and A distributions includes the finite regular terms not included in A and 
logarithmic terms from the final-state fragmentation singularities, with the latter subtracted 
when Qt < £™, as described in Sec. Ill CI The data clearly disfavor the fixed-order prediction 
in the region of low Qt- 

Figure Mb) features the resummed W + Y contribution (solid curve). Resummation of 
the initial-state logarithmic terms renders W finite in the region of small Qt- The sum of W 
and Y includes the resummed initial-state singular contributions plus the remaining relevant 
terms in P. Since P provides a reliable fixed-order estimate at large Qt, we present our 
final resummed prediction by switching from W + Y to P at the point at which the two 
differential cross sections (as functions of Q, Qt and y) cross each other. In contrast to the 
fixed-order (dashed) curve P in Fig. E](b), the agreement with data is improved at the lowest 
values of Qt, where resummation brings the rate down, and for Qt = 12 — 32 GeV, where 
the resummed logarithmic terms increase the rate. 

The resummed predictions for the Tevatron experiments are practically insensitive to the 
choice of the resummation scheme and the nonperturbative model 0. About 75% (25%) of 
the total rate at the Tevatron with CDF cuts imposed comes from the qq + qg + qg (gg + gqs) 
initial state. The fractions for the cuts used by D0 are 84% and 16%. The gg + gqs 
contribution falls steeply after Qt > 22 GeV, because the gluon PDF decreases rapidly with 
parton fractional momentum x [3j. 

The distribution in the difference A(p of the azimuthal angles of the photons is shown 
in Fig. El As is true for the transverse momentum distribution in the limit Qt — > 0, the 
distribution computed at fixed order is ill-defined at Atp = ir. The resummed distribution 
shows a larger cross section near Atp = 2.5 rad, in better agreement with the data. In the 
region of small Atp < 7r/2, the fixed-order and the resummed predictions are the same, a 
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Figure 6: The difference in the azimuthal angles of the two photons in the laboratory frame 
predicted by the resummed (solid) and fixed-order (dashed) calculations, compared to the CDF 
data. 



result of our matching of the resummed and fixed-order distributions at mid to high values 
of Qt- Although the cross section is not large in the region A<p < ir/2, there is an indication 
of a difference between our predictions and data in this region, a topic we address below. 



3. The region Qt > Q 



It is evident from Fig. [3] that the Atp < n/2 region is populated mostly by events with 
Qt > Q- New types of radiative contributions may be present in this region, including various 
fragmentation contributions described in Sec. Ill CI and enhancements at large | cos 0*1 in the 
direct production rate. 

While experimental isolation generally suppresses long-distance fragmentation, a greater 
fraction of fragmentation photons are expected to survive isolation when A<p < n/2. Besides 
single-photon 'one-fragmentation' and 'two-fragmentation' contributions (with one photon 
per fragmenting parton), one encounters additional logarithmic singularities of the form 
\og{Q / Qt) ■ We noted in Sec. Ill Cl that these logarithms are associated with the fragmentation 
of a parton carrying large transverse momentum Qt into a system of small invariant mass Q 
28l |. a light 77 pair in our case. Small-Q 77 fragmentation of this kind is not implemented 



yet in theoretical models. Therefore, we are prepared for the possibility that both the 
fixed-order calculation and our resummed calculation may be deficient in the region Qt ^> 
Q. A detailed experimental study of the region Qt > Q may offer the opportunity to 
measure the parton to two-photon fragmentation function D 11 {z\^ Z2), provided that the 
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single-photon 'one- fragmentation' function D y (z) is determined by single-photon data, and 
the low-Q logarithmic terms are properly resummed theoretically. 

In addition to the \ow-Q fragmentation, the small-Ay? region may be sensitive to large 
higher-order contributions associated with t- or u-channel exchanges in the qq — > 77 and 
99 77 subprocesses. In the Born processes in Figs. [H(a) and (h), the t- and w-channel 
singularities arise at cos6>* « ±1 and A(p « n. These singularities are excluded by the p T 
cuts in Eq. (1231 ). but related residual enhancements in the NLO contributions may still persist 
at I cos^l ~ 1 and A(p — > 0, not excluded by the cuts (cf. Fig. [3]) . Because | cos 6**1 is large in 
such events, they tend to have substantial \Ay\, so they are retained by the Ai? 77 > 0.3 cut. 
In contrast, the low-Q fragmentation contributions tend to be abundant at small \Ay\. It 
may be therefore possible to distinguish between the large- 1 cos#*| and fragmentation events 
at small A(p based on the distribution in \Ay\. 

We expect much better agreement of our predictions with data if the selection Qt < Q is 
made. This selection preserves the bulk of the cross section and assures that a fair comparison 
is made in the region of phase space where the predictions are most valid. 



4- Fragmentation and comparison with the DIPHOX code 

One way to obtain an estimate of theoretical uncertainty is to compare theoretical ap- 
proaches in various parts of phase space, including small A(p. We handle the collinear 
final-state photon singularities in the manner described in Sec. UH without including photon 
fragmentation functions explicitly. An alternative calculation implemented in the DIPHOX 
code 3] includes NLO cross sections for single-photon fragmentation processes. Neither 



code includes a term in which both photons are fragmentation products of the same final- 
state parton, i.e., the diphoton fragmentation function D 11 (zi, Z2). 

In Ref. [1] we show comparisons of our predictions with those of DIPHOX along with the 
CDF data. Here in Fig. [TJ we show analogous plots of the invariant mass and transverse mo- 
mentum distributions for D0 cuts. We note that our fixed-order qq + qg contribution agrees 
well with the direct contribution in DIPHOX. This agreement is particularly impressive in 
the region of large Qt, where both codes use the same fixed-order formalism to handle direct 
contributions. A contribution from the gg channel is also present in both codes, computed 
at LO in DIPHOX but at NLO+NNLL in our case. Since the gg + gqs contribution is 
not dominant (especially in the high Qt region), this difference does not have a significant 
impact on the comparison. 

The explicit single-photon fragmentation contributions in DIPHOX (mostly 'one- 
fragmentation' contribution) are quite small for the nominal hadronic energy ~ 1 GeV 
in the cone around each photon. Exceptions occur in the region Qt < E^f°, where the 
fragmentation contributions to da/dQT have logarithmic singularities, and in the Aip — > 
region, where fragmentation is comparable to the direct contributions. Our isolation pre- 
scription reproduces the integrated DIPHOX rate well for < Qt < E^°, leading to close 
agreement between the resummed and DIPHOX inclusive rates for most Q values. 

Returning to the CDF measurement, we remark that the resummed and DIPHOX cross 
sections for the same Elf° = 1 GeV underestimate the data within two standard deviations 
for Q < 27 GeV, Q T > 25 GeV, and Aip < 1 rad (cf. the relevant figures in Ref. Q). 
The DIPHOX cross section can be raised to agree with data in this "shoulder" region, if 
a much larger isolation energy {E™ = 4 GeV) is chosen, and smaller factorization and 
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Figure 7: Comparison of our resummed and DIPHOX predictions for (a) the invariant mass and 
(b) transverse momentum distributions of 77 pairs for D0 kinematic cuts. The solid curves show 
our resummed distributions with all channels included. The dashed and dotted curves illustrate the 
resummed and DIPHOX distributions in the qq + qg channel. 



renormalization scales are used (/xp = [Ir = Q/2). These are the choices made in the CDF 
study [H]. Since is an approximate characteristic of the experimental isolation, one 
might argue that both Etf° = 1 and 4 GeV can be appropriate in a calculation to match 
the conditions of the CDF measurement. The direct contribution is weakly sensitive to E™°, 
while the one-fragmentation part of da/ 'dQ T is roughly proportional to E™° (cf. Section HTCl) . 
The one-fragmentation contribution is enhanced on average by 400% if E™ is increased in 
the calculation from 1 to 4 GeV. The rate in the shoulder region is enhanced further if the 
factorization scale \ip is reduced. 

Since the theoretical specifications for isolation and for the fragmentation contribution 
are admittedly approximate, we question whether great importance should be placed on the 
agreement of theory and experiment in the region of small Atp or in the shoulder region 
in the Qt distribution. A straightforward way to reduce sensitivity to fragmentation is to 
require Q > 27 GeV or Qt < Q, as discussed above. The two cuts have similar effects on 
the event distributions. Figure [8] shows the effects of the Qt < Q cut on the Qt and Atp 
distributions. The cut Qt < Q is particularly efficient at suppressing the fragmentation Qt 
shoulder (and the region of small Aip altogether), while only a small portion of the event 
sample is lost. This cut is especially favorable, since it constrains the comparison with data 
to a region where the theory is well understood and has a small uncertainty. Furthermore, 
with the requirement of Qt < Q, the dependence of differential cross sections on the choices 
of isolation energy E^° and factorization scale \ip is greatly reduced to the typical size of 
higher-order corrections. We predict that if a Qt < Q cut, or a Q > 27 GeV cut, is applied 
to the Tevatron data, the enhancement at low Aip and intermediate Qt associated with the 
fragmentation contribution will disappear. This is an important conclusion of our study, and 
we urge the CDF and D0 collaborations to apply these cuts in their future analyses of the 
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Figure 8: Predicted cross sections for diphoton production in pp — ► 77 V at y S = 1.96 TeV as a 
function of (a) the 77 pair transverse momentum Qt and (b) the difference A(p in the azimuthal 
angles of the two photons. Our resummed predictions (solid) are shown together with DIPHOX 
predictions for the default isolation energy E^f° = 1 GeV and factorization scale fiF = Q (dashed), 
and for = 4 GeV, [if = Q/2 (dotted). We impose the condition Qt < Q to reduce theoretical 
uncertainties associated with fragmentation. 



diphoton data. 



5. Average transverse momentum 

An important prediction of the resummation formalism is the change of the transverse 
momentum distribution with the diphoton invariant mass. This dependence comes, in part, 
from the hiQ 2 dependence in the Sudakov exponent, Eq. (fl"7l) . and it is desirable to identify 
this feature amid other influences. In Fig. [9](a), we show normalized resummed transverse 
momentum distributions for various selections of the invariant mass of the photon pairs. 
Without kinematical constraints on the decay photons, the Qt distribution is expected to 
broaden with increasing Q, and the position of the peak in da/dQT to shift to larger Qt 
values. The shift of the peak may or may not be observed in the data depending on the 
chosen lower cuts on pt of the photons, which suppress the event rate at low Q and Qt- The 
interplay of the Sudakov broadening of the Qt distribution and kinematical suppression by 
the photon pt cuts is reflected in the shape of da/dQT in different Q bins. 

According to dimensional analysis, the average {Qt) in the interval Qt < Q may be 
expected to behave as 

{Qt)q t <q = Qf(Q/VS), (25) 
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Figure 9: (a) Resummed transverse momentum distributions of photon pairs in various invariant 
mass bins used in the CDF measurement, normalized to the total cross section in each Q bin. (b) 
The average Qt as a function of the 77 invariant mass, computed for Qt < Q- 



where the scaling function f(Q/y/S) reflects phase space constraints, dependence on the 
Sudakov logarithm, and the x dependence of the PDFs. Figure [9](b) shows our calculated 
diphoton mass dependence of (Qt)- The linear increase shown in Eq. (l25|) is observed over 
the range 30 < Q < 80 GeV. For values of Q below the kinematic cutoff at about 30 GeV, the 
cuts shown in Fig. 3 suppress diphoton production at small Qt, and (Qt/Q) grows toward 
1 as Q decreases (corresponding to production only at Qt close to Q). For Q ~ 80 GeV 
and above, we see a saturation of the growth of (Qt), a reflection of the influences of the 
x dependence of the PDFs and other factors. Similar saturation behavior is observed in 



calculations of (Qt) in other processes [36J]. It would be interesting to see a comparison of 



our prediction with data from the CDF and D0 collaborations. 



B. Results for the LHC 



1. Event selection 



To obtain predictions for pp collisions at the LHC at v S = 14 TeV, we employ the cuts 
on the individual photons used by the ATLAS collaboration in their simulations of Higgs 



boson decay, h — > 77 [371] ■ We require 



transverse momentum p T > 40 (25) GeV for the harder (softer) photon, 
and rapidity |y 7 | < 2.5 for each photon. 



(26) 
(27) 



In accord with ATLAS specifications, we impose a looser isolation restriction than for our 
Tevatron study, requiring less than = 15 GeV of hadronic and extra electromagnetic 
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Figure 10: Resummed transverse momentum distributions of photon pairs in various invariant 
mass bins at the LHC. The cuts listed in Eqs. ([26]) and (|27l) are imposed. The Qt distribution for 
70 < Q < 115 GeV with an additional constraint Qt < Q is shown as a dotted line. 



transverse energy inside a AR = 0.4 cone around each photon. We also require the separation 
Ai? 77 between the two isolated photons to be above 0.4. 

The cuts listed above, optimized for the Higgs boson search, may require adjustments in 
order to test perturbative QCD predictions in the full 77 invariant mass range accessible at 
the LHC. The values of the p T cuts on the photons in Eq. ([26]) preserve a large fraction of 
Higgs boson events with Q > 115 GeV. These cuts may be too restrictive in studies of 77 
production at smaller Q, considering that the two final-state photons most likely originate 
from a 77 pair with small Qt and have similar values of p T of about Q/2. The pt cuts 
interfere with the expected Sudakov broadening of Qt distributions with increasing diphoton 
invariant mass, as discussed in Section Ull A 51 We further note that the asymmetry between 
the p T cuts on the harder and softer photons is necessary in a fixed-order perturbative QCD 
calculation, but it is not required in the resummed calculation. At a fixed order of a s , 
asymmetry in the p T cuts prevents instabilities in da/dQ caused by logarithmic divergences 
in da/dQT at small Qt- Such instabilities are eliminated altogether once the small-Q^ 
logarithmic terms are resummed to all orders of a s . Here we do not consider alternative 
p\ cuts, although experimental collaborations are encouraged to employ relaxed and/or 
symmetric cuts to increase the 77 event sample in their data analysis. 



2. Resummed Qt distributions and average transverse momentum 

Figure [10] shows transverse momentum distributions of the photon pairs for various in- 
variant masses. The average 77 transverse momentum grows with Q, as demonstrated by 
Fig. [TH However, the rate of the growth decreases monotonically with Q, for similar reasons 
as at the Tevatron. 

The 77 distributions in Q and A<p for different combinations of scattering subchannels 
and choices of theoretical parameters are discussed in Refs. [1, 0|. In all ranges of Q, the 77 
production rate is dominated by a large qg contribution, accounting for about 50% of the 
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Figure 11: The average Qt at the LHC as a function of the 77 invariant mass Q. 



fixed-order (NLO) rate. Although this number depends on the choice of the factorization 
scheme and scale, and, on the other hand, separate treatment of the qq and qg cross sections 
is not meaningful in the resummation calculation (§], it nonetheless reflects, in a crude way, 
the increased relative importance of the qg cross section. The gg + gqs channel contributes 
about 25% at Q ~ 80 GeV (the location of the cutoff in da/dQ due to the cuts on p T ) and 
less at larger Q. As at the Tevatron, the dependence of the cross sections on the resummation 
scheme is small (jjj. The dependence on the nonperturbative model can also be neglected, as 
long as the nonperturbative function does not vary strongly with x 0|. 



3. Final-state fragmentation and comparison with DIPHOX 

The impact of the final-state fragmentation at the LHC can be evaluated if we compare 
our results with DIPHOX predictions. The transverse momentum and invariant mass dis- 
tributions in the qq + qg channel from the two approaches are shown in Fig. [T2l In both 
calculations, quasi-experimental isolation removes direct NLO events with collinear final- 
state photons and partons when Qt > E™ = 15 GeV, but not when Qt is below E^°. 

Concentrating first on 77 events with Qt > E^°, we observe that, at Qt > 80 GeV, the 
resummed qq+qg cross section reduces to the direct fixed-order cross section, evaluated in the 
same way as in the DIPHOX code. Our resummed and the direct DIPHOX cross sections, 
shown as solid and dashed curves, respectively, in Fig. [121(a) consequently agree well at 
large Qt- At smaller Qt, the resummed cross section is enhanced by towers of higher-order 
logarithmic contributions. On the other hand, the full qq + qg DIPHOX rate (shown as a 
dotted line) also includes single-photon fragmentation contributions, which add to the direct 
production cross section. For the nominal isolation parameters, the explicit fragmentation 
contribution constitutes about 25% of the full DIPHOX rate for 60 < Q T < 120 GeV. Its 
magnitude increases approximately linearly with the assumed E^° value. 

For Qt < E^°, the final-state collinear region of the direct contribution is regulated 
by the collinear subtraction prescription adopted in the resummation calculation, whereas 
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Figure 12: Transverse momentum and invariant mass distributions da/dQ in the qq + qg channel 
obtained in the resummation (solid) and DIPHOX (dotted) calculations. 



the fragmentation singularity is subtracted from the direct contribution and replaced by 
photon fragmentation functions in the DIPHOX calculation. Subtraction of singularities in 
DIPHOX introduces integrable singularities in da/dQx at different values of Qt below E™°- 
The origin of the final-state logarithmic singularities at values of Qt below E l ^° is discussed 
in Refs. 22, 23, 24j. For Qt < E^°, the DIPHOX curves represent the average over singular 



contributions in this Qt interval. 

After the fragmentation singularity is subtracted, the DIPHOX direct contribution 
(dashed line) is on average below our resummed qq + qg rate (solid line) over most of the 
range of Qt shown in Fig.[T2Ta). After integration over all Qt, our resummed and DIPHOX 
qq + qg cross sections agree within 10-20% at most values of Q (cf. Fig. [121(b)). with our 
resummed rate being below the DIPHOX rate at all Q. The largest difference occurs at the 
lowest values of Q (below the cutoff), where the rates can differ by a factor of 2. In this 
region, corresponding to diphoton events with small A(p and Qt larger than Q, the photon 
fragmentation contributions included in the DIPHOX calculation are large in comparison to 
the direct rate. Finally, we note that the integrated rate in DIPHOX is more stable with 
respect to variations in than the differential distributions in DIPHOX, b ecause Erp° 
dependence for Qt > Elf° is canceled to a good degree by E^° dependence for Qt < E^°. 

To obtain the final 77 production cross sections, after inclusion of all channels, we combine 
the respective qq + qg results with the resummed NLO gg + gqs cross section in our case and 
with the LO gg cross section in the DIPHOX case. The distributions in the 77 invariant 
mass Q, the transverse momentum Qt, and the azimuthal angle separation Aip in the lab 
frame are shown in Fig. [T"3l For the cuts chosen, the LO gg and the resummed gg + gqs total 
rates constitute about 9% and 20% of the total rate. The resummed and DIPHOX invariant 
mass distributions (Fig. [131(a)) are brought closer to one another as a result of the inclusion 
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Figure 13: Invariant mass, transverse momentum, and Aip distributions from our resummed cal- 
culation and from DIPHOX at the LHC. We show our fixed-order (dashed) and resummed (solid) 
distributions. All initial states are included in both calculations, and the single-7 fragmentation 
contributions are included in DIPHOX. 
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Figure 14: Invariant mass and transverse momentum distributions from our resummed, NLO, and 
DIPHOX calculations at the LHC, with the Qt < Q constraint imposed. 

of the gg + gqs contribution in the resummed calculation. For Qt ^ 0, the full DIPHOX Qt 
distribution in Fig. [131(b) is determined entirely by direct plus fragmentation contributions 
(the same as in Fig. [121(a)). because the LO gg cross section contributes at Qt = only. In 
contrast, our resummed gg + gqs contribution modifies the event rate at all Qt- 

The resummed and DIPHOX rates are in a reasonable agreement for 1.5 < A<p < 2.5, as 
shown in Fig. flBTc). In the A(p — > ti limit, the fixed-order rates in DIPHOX diverge because 
of the singularities at small Qt, while our resummed rate yields a finite value. For A(p < 1.5, 
the DIPHOX cross section is enhanced by photon fragmentation contributions. As at the 
energy of the Tevatron, theoretical uncertainties are greater at small A<p. 

Predictions are most reliable when Qt < Q (and the angles 9* and <p* are away from 
or it). With the Qt < Q cut imposed, the uncertain large-Qy photon fragmentation 
contributions are suppressed, and the resummed and DIPHOX cross sections agree well at 
large Q T (cf. Fig. (11(b)). The Q T distribution in the interval 70 < Q < 115 GeV with the 
Qt < Q constraint is shown in Fig. [TO] by a dotted curve. Distributions in the other two 
mass bins in Fig. [10] are essentially not affected by this cut in the Qt range presented. 

Our calculation captures the dominant contributions to 77 production at the LHC. How- 
ever, as we noted, direct qg scattering, evaluated at order 0(a s ) in our calculation, is the 
leading scattering channel in the region relevant for the Higgs boson search at the LHC. 
It is important to emphasize that the final-state collinear radiation is not the main reason 
behind the enhancement of the qg rate, which is increased predominantly by contributions 
from non-singular phase space regions. Consequently, the qq + qg direct rate is only weakly 
sensitive to adjustments in the isolation parameters EJf° and AR [l(|. The unknown 0{a 2 s ) 
contributions to qg scattering may be non-negligible, and it would be valuable to compute 
them in the future when LHC data are available. 
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Figure 15: Comparison of the normalized Higgs boson signal and diphoton background distributions 
at the LHC, both computed at NNLL accuracy. The Higgs boson mass is taken to be mg = 130 
GeV, and the background is calculated for 128 < Q < 132 GeV. 

C. Comparison with Higgs boson signal distributions 

We highlight some similarities and differences between the production spectra for the 
Higgs boson signal and the QCD background discussed in this paper. We focus on the 
diphoton decay mode of a SM Higgs boson produced from the dominant gluon-fusion mech- 
anism, gg — > h° —>■ 77, where the Higgs boson production cross section is calculated at the 
same order of precision as the QCD continuum background. We include initial-state QCD 
contributions at 0(af) (NLO) and resummed contributions at NNLL accuracy. These con- 



tributions are also coded in ResBos [38|], and we can apply the same cuts on the momenta 
of the photons to the signal and background. Our findin gs s hould remain broadly applicable 
after the NNLO corrections to Higgs boson production 139, S3 are included. We compute 



the background in the range 128 < Q < 132 GeV, and the signal at a fixed Higgs boson 
mass m H = 130 GeV. We impose the kinematic selection Qt < Q-, but its influence is not 
important at the large values of diphoton mass of interest here. 

The cross section times branching ratio for the Higgs boson signal is substantially smaller 
than the QCD continuum. To better illustrate their differences, Fig. [15] presents distributions 
normalized to the respective total rates. The top-left panel shows normalized transverse 
momentum distributions of photon pairs. The signal and background peak at about 12 and 
5 GeV, respectively. The average values of Qt are 26 and 23 GeV, computed over the range 
to 75 GeV. 

Differences in the shapes of these Qt spectra can be attributed to the different structure 
of the leading terms in the initial-state Sudakov exponents and to the effects of final-state 
photon fragmentation. The Higgs boson signal is controlled by the characteristics of the 
gg+gqs initial state, whereas the continuum is controlled primarily by the qq+qg initial state. 
Because the dominant Sudakov coefficient Aq oc Cp in the qq case is smaller than A g ® = 
(C A /C F )A ( q ) in the gg case, the resummed qq + qg initial-state radiation produces narrower 
Qt distributions than gg + gqs initial-state radiation. About 80% of the diphoton rate is 
provided by the qq + qg channel, implying a narrower Qt distribution of the background, if 
based on the value of A^ alone. 

The continuum background contribution is also enhanced by final-state radiation in qg 
scattering. The Qt profile of the final-state collinear terms depends more on the isolation 
model (including and AR) than on the initial-state Sudakov exponent. For the nominal 
ATLAS cuts, the final-state collinear contribution in our calculation hardens the background 
Qt distribution, diminishing its difference from the Higgs boson signal distribution. More 
effective isolation may reduce the impact of the final-state radiation on Q T distributions. 

Another difference between the signal and continuum is observed in the distribution in 
the azimuthal angle of the photons, such as the angle ip* in the Collins-Soper frame shown 
in the top-right panel of Fig. [151 This distribution is qualitatively the same if integrated 
over all Qt, as in Fig. [IH or integrated above some minimal Qt value, as in an experimental 
measurement. Without isolation imposed, the spin-0 Higgs boson signal must be flat in </?*, 
but the QCD background peaks toward ip* = and re (i.e., sin^* = 0) as a result of the 
final-state qg singularity. 2, 3 Isolation cuts suppress both the signal and the background for 
sin 92* < sin AR. The result is a signal distribution with a broad peak near = ir/2, while 
the background favors values of near and ir. A selection of events with in the vicinity 
of 7r/2, and Qt large enough, helps to reduce the impact of the qg background. In the lab 
frame, a related distribution is in the variable ~ <^4t| , where ipiT is the azimuthal angle 
between pj? and Qt- The signal (background) processes tend to have more events with large 



2 By definition, the recoil parton 5 always lies in the Oxz plane (has zero azimuthal angle) in the Collins- 
Soper frame. For the final-state singularity to occur at NLO, the photons should be in the same plane 
with 5, i.e., have sin tp* = 0. 

3 One of the resummed structure functions for the gg background is modulated by cos2(^» (see Sec. Ill B| . 
but we neglect this modulation in our present calculation. 



(small) magnitude of 1^3^ — <p^r\. 

A third potential discriminator between the signal and background is the difference in 
the rapidities Ay = yhard — Vsoft of the photons with harder and softer values of p T m the 
lab frame, calculated on an event by event basis. This distribution is displayed in the lower- 
left frame of Fig. [15j The background distribution peaks at the origin, while the signal is 
almost flat over a wide range of Ay. Different spin correlations in the decay of a spin-0 
Higgs boson from those characteristic of QCD background processes are the source of this 
distinction. Discrimination based on this difference can improve the statistical significance 
of the signal We note that our resummed calculation does not exhibit the kinematic 
singularity at Ay « 2 present in the finite-order cross section and obvious in Fig. 10 of 
Ref. [10], where the distribution with respect to y* = Ay/2 is shown. The discontinuity in 
da /dQx caused by the finite-order approximation is resummed in our calculation, yielding a 
smooth result. 

The rapidity difference is related to the scattering angle in the Collins-Soper frame: 
tanh(Ay/2) = cos9* when Qt is zero. The cos9* distribution is shown in the lower-right 
frame of Fig. [T51 The difference between the signal and background rates is even more 
pronounced in this variable, clearly expressing the difference in the spin correlations of the 
systems producing the photons. 

A comparison of Qt distributions in the top-left panel of Fig. [15] suggests that the signal 
versus background ratio would be enhanced if a cut is made to restrict Qt > 10 GeV. 
After applying this cut, we may again examine the distributions in the rapidity difference of 
the two photons, the scattering angle in the Collins-Soper frame, and the azimuthal angle 
distribution of the photons in the Collins-Soper frame. The results are qualitatively similar 
to those in Fig. [Tj3 and are not shown here. A more efficient procedure to increase the 
Higgs boson discovery significance is to apply a simultaneous likelihood analysis to several 
kinematic distributions. Based on the present discussion, we would argue that the resummed 
Qt, v 9 *! an d cos#* distributions are good discriminators between the Higgs boson signal and 
background in such an analysis. 



IV. CONCLUSIONS 

The theoretical study of continuum diphoton production in hadron collisions is interesting 
and valuable for several reasons: there are data from the CDF and D0 collaborations at 
Fermilab with the promise of larger event samples; there are new theoretical challenges 
associated with all-orders soft-gluon resummation of two-loop amplitudes; and continuum 
diphotons are a large standard-model background above which one may observe the products 
of Higgs boson decay into a pair of photons at the LHC. 

In this paper and Refs. [3, S], we present our calculation of the fully differential cross 
section da / (dQdQTdydQ*) as a function of the mass Q, transverse momentum Qt, and 
rapidity y of the diphoton system, and of the polar and azimuthal angles of the individual 
photons in the diphoton rest frame. Our basic QCD hard-scattering subprocesses are all 
computed at next-to-leading order (NLO) in the strong coupling strength a s , and we include 
the state-of-art resummation of initial-state gluon radiation to all orders in a s , valid to next- 
to-next-to-leading logarithmic accuracy (NNLL). Resummation is essential for a realistic and 
reliable calculation of the Qt dependence in the region of small and intermediate values of 
Qt, where the cross section is greatest. It is also needed for stable estimates of the effects of 



30 



experimental acceptance on distributions in the diphoton invariant mass and other variables 



Our analytical results are included in a fully updated ResBos code [31|, |32fl. This nu- 
merical program allows us to impose selections on the transverse momenta and angles of the 
final photons, in order to match those employed by the CDF and D0 collaborations, as well 
as those anticipated in experiments at the LHC. Our predictions are especially pertinent in 
the region Qt < Q- We show that our results at the Tevatron and at the LHC are insensitive 
to the choice of the resummation scheme and of the nonperturbative functions required by 
the integration into the region of large impact parameter. 

The published collider data are presented in the form of singly differential distributions. 
We follow suit in order to make comparisons, and we find excellent agreement with data, as 
shown in Sec. IHH We recommend that more differential studies be made, and, to motivate 
these, we present predictions for the changes expected in the Qt distribution as a function 
of mass Q, and for the dependence of the mean transverse momentum on Q. 

We make predictions for continuum diphoton mass, transverse momentum, and angular 
distributions at the energy of the LHC. Moreover, we contrast in Fig. [15] the shapes of some 
of these distributions with those expected from the decay of a Higgs boson. The distinct 
features of the signal and background suggest that that the Higgs boson discovery significance 
can be increased via a simultaneous likelihood analysis of several kinematic distributions, 
particularly the resummed Qt, y?*, and cos#* distributions. 

Another approach to the computation of continuum diphoton production is presented 



by the DIPHOX collaboration This calculation includes both the direct production 

of photons from hard-scattering processes and the photons produced from fragmentation of 
(anti-)quarks or gluons. It is valid at NLO, except for the gg subprocess, which is included 
at leading order only. The DIPHOX code is useful for rates integrated over transverse mo- 
mentum, but it is not designed to predict the transverse momentum distribution or other 
distributions sensitive to the region in which the transverse momentum of the diphoton pair 
is small. Compared to a fixed-order calculation, such as direct photon pair calculation in 
DIPHOX, our calculation improves the theoretical prediction for event distributions which 
are sensitive to the region of low Qt- Furthermore, our calculation includes the NLO con- 
tribution from the combined gg + gqs channel, leading to more accurate predictions at the 
LHC, where the gg + gqs contribution is generally not small. 

Only isolated, not inclusive, photons are identified experimentally. While it is straightfor- 
ward to define an isolated photon in a given experiment, it is challenging to devise a theoret- 
ical prescription that can match the experimental definition, short of first understanding the 
long-distance dynamics of QCD. The isolated diphoton production rate is modeled in the 
DIPHOX code by including explicit photon fragmentation function contributions at NLO 
accuracy. A shortcoming of this approach (as well as of our method for treating isolation) 
is that one cannot accurately represent photon fragmentation without including final-state 
parton showering in the presence of isolation constraints. There is inevitable ambiguity and 
uncertainty in the choice of the "isolation energy" used to define an isolated photon the- 
oretically for comparison with the isolated photon measured experimentally. As shown in 
Sec. IIIIl the DIPHOX cross section can vary by a large factor in some regions of phase space 
at the Tevatron when is changed from 1 GeV to 4 GeV. 

Our approach is to concentrate on physical observables which are less sensitive to the frag- 
mentation contributions. We apply the "collinear subtraction" prescription or the "smooth- 
cone isolation" prescription to define an isolated photon in our calculation. We find good 
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agreement with the data, except in the region with small Q and Aip < n/2, consistent with 
our theoretical expectation that higher-order direct photon production and photon fragmen- 
tation contributions can strongly modify the rate of diphoton pairs in this region. We suggest 
that much better agreement with current and future data will be obtained if an addition 
requirement of Qt < Q is applied. With this cut, the fragmentation contributions are largely 
suppressed. With the cut Qt < Q cut applied to the Tevatron data, the enhancement at 
low Aip and intermediate Qt (the shoulder region) should disappear. We urge the CDF and 
D0 collaborations to apply these cuts in future analyses of the diphoton data. 

In our calculation, we identify an interesting spin-flip contribution (with cos2<£>* depen- 
dence) in the gg channel, cf. Ref. [3J, and we suggest that measurements be made of the 
distribution of (p* as a function of Qt- All-orders resummation of the gluon spin-flip contri- 
bution may be needed when a larger statistical sample of diphoton data is available. 

The contributions from qg + qg processes become more important at the LHC than at the 
Tevatron, and calculations at a higher order of precision may be warranted eventually. To 
improve the theoretical prediction in the region of phase space with Qt < and ~ 
or 7r, a joint resummation calculation is needed in which the effects of both the initial- and 
final-state multiple parton emissions are treated simultaneously. 

Although we emphasize that better agreement of our predictions with data should be 
apparent if the selection Qt < Q is made, we also point out that the region Qt > Q should 
manifest very interesting physics of a different sort. Additional logarithmic singularities of 
the form \og(Q / Qt) are encountered in the region Qt Q. These logarithms are associated 
with the fragmentation of a parton carrying large transverse momentum Qt into a system 
of small invariant mass Q [23, 2§|, a light 77 pair in our case. Small-Q 77 fragmentation 



of this kind is not implemented yet in theoretical models. Experimental study of the region 
Qt ^ Q may offer the opportunity to measure the parton to two-photon fragmentation 
function D 71 (zi, Z2). 
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Appendix A: SUMMARY OF PERTURB ATIVE COEFFICIENTS 

In this appendix we present an overview of the perturbative QCD expressions for the 
resummed and asymptotic cross sections used in our computation. 

The functions A a (Ci,p), B a (C±, C2, p), C a / ai {x, b; C1/C2, n), and h a (Q,9*) are introduced 
in Sec. OH These functions are derived as perturbative expansions in the small parameter 
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71=1 
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C a/ai U &5 ^, A«) = £ Ci)l (x, bp, £■) (^) " ; K(Q, 6.) = (9, 

^ 2 / n= Q 2 \ / n=0 



«.(Q)\ n 



7T 



The value of a perturbative coefficient for a set of scales C\jb and C 2 Q can be expressed 
in terms of its value F^ n ' c ^ obtained for the "canonical" combination C\ = Cq and C 2 = 1. 
Here c = 2e~ 7E ~ 1.123, where 7# = 0.5772... is the Euler constant. The relationships 
between and F^ n ' c ' take the form 



^ 1} (Ci) = 

4 2) (^i) = ^-4^0 In 00 

4 3) (Ci) = 4 3 ' c) "24 2 ' c) /5oln 
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*Ja | '-'0 2 



(l>c) / „2^2 
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They depend on the QCD beta-function coefficients fa = (11N C - 2N f )/6, A = (17A 2 
5A C A/ - 3C F N f )/6 for A c colors and N f active quark flavors, with C F = (A 2 - 1)/(2A C ) 
4/3 for A c = 3. The relevant 0(a s ) splitting functions P a / ai (x) are 



p q / q 


= c F ( 


l + z 2 \ 
l-x 


' Pq/g 


p 9 / 9 


= 2C A 
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1 — X 
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[1 + 2X + 2X 2 ); P 9/gs =C F 



[1-xY + l 



x 



+ x(l — x) 



+ /3o6(l-x). 



(A7) 
(A8) 



The coefficients hP)(9*), B^ 2 \ and & x > depend on the resummation scheme. The hard- 
scattering function is 



« s (g)v a (^) 



7T 



4 



(A9) 



where <5 S = in the CSS scheme and <5 S = 1 in the CFG scheme. The functions V q {9 tt ) for 
qq —> 77 scattering and V g (0*) for gg — > 77 scattering are derived in Refs. 12] and [13], 
respectively. 
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For the qq + qg initial state, we obtain the following expressions for the coefficients A, B, 
and C: 



4 1, c) = CV; 

4 2, c) = c F 



Y67 


7T 2 > 










V36 ~ 


12, 



(A10) 



Af c) 



3 



H(2,c) 



C^iV/ 55 



3 



/ 108 A V 24 7 20 216 
7C(3) 5tt 2 209 



245 
"96" 



Cf', 



C F 2 



12 



3 _ 7T 2 

8 2 
1 2tt 2 



108 432 



6C(3) 



( — + ^-3C(3) 
1 24 18 SV ; 



(i,c)/__% I Cf 



9 

(°), 



C F N f T R ( - 
C$(s) = 8 jk 8(l-x); q%(x) = 

1 



C F 7T 2 

12 
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C f/g\ X ) = ^(1-X) 



(All) 



Here Ca = N c , Tr = 1/2, and the Riemann constant C(3) = 1.202 .... The C functions are 
given for j, k = u,u,d,d, . . .. These coefficients are taken from EES]. 
Similarly, the A, B, and C coefficients in the gg + gqs channel are 



A k ' c) = (C A /C F )A^ C \ for k — 1,2, 3; 
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These coefficients are taken from Refs. [12j, LL3J, |44|, 145 



Appendix B: COMPONENTS OF THE ASYMPTOTIC CROSS SECTIONS 



In Sec. Ill Bl we introduce asymptotic small-Qr approximations for the qq + qg and gg + gqs 
NLO cross sections, 



Aqq(Q, Q T , y, &*) 
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i=u,u,d,d,... 
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(Bl) 
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and 



A gg (Q,Q T ,y,n*) = -<jE g (0*) 5(Q T )F gtS (Q,y,9,) + F gi+ (Q,y,Q T ) 



+Z g (e*, ( p*)F' g (Q,y,Q T 
The functions F in these equations are defined as 

Fi,s{Q,y,d*) = f q . l /h 1 { x i^F)U/h 2 {x2,^F) (l + 2—h£ ) (e*) 
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Expressions for the coefficients .Ai , B a , ha (9*), C^/(x), and splitting functions P a / C (x), 
are listed in Appendix lAl Summation over all relevant parton flavors a' = g,u,u,d,d, ... 
for a = q and a' = g, qs for a = g is assumed. In addition, the </vdependent part 
£'„(#*, (0*)F'(Q, y, Qt) of the gg + gqs asymptotic cross section A gg contains a splitting 



function 



P; g (x) = 2C A (1 - x)/x, 



(B8) 



contributed by the interference of splitting amplitudes with opposite gluon polarizations in 
the helicity amplitude formalism [46j, l47|, |48|, [49||. The origin and behavior of this spin-flip 
function are discussed in Ref. (§]. 
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